Temperature-Polarization Correlations from Tensor Fluctuations 
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We study the polarization-temperature correlations on the cosmic microwave sky resulting from an 
initial scale invariant spectrum of tensor (gravity wave) fluctuations, such as those which might arise 
during inflation. The correlation function has the opposite sign to that for scalar fluctuations on 
large scales, raising the possibility of a direct determination of whether the microwave anisotropics 
have a significant tensor component. We briefly discuss the important problem of estimating the 
expected foreground contamination. 



COBE's detection of microwave background tempera- 
ture anisotropies W has excited great interest in theories 
of cosmic structure formation. Theoretical mechanisms 
for producing CMB anisotropies include primordial en- 
ergy density (scalar) and gravitational wave (tensor) fluc- 
tuations generated during inflation, or the gravitational 
flelds induced by cosmic defects. These three mechanisms 
make surprisingly similar predictions for anisotropies on 
the large angular scales probed by COBE, and further 
observations are needed in order to discriminate between 
them. It is clearly of particular value to identify dis- 
tinctive signals associated with specific physical effects, 
in order to determine in as direct a manner as possible 
which mechanisms produced these anisotropies. 

The level of temperature fluctuations on smaller scales 
provides the simplest such test. For example, within 
inflationary models, the anisotropies due to adiabatic 
scalar fluctuations generally increase on scales of or- 
der the horizon at last scattering (the "Doppler peak") 
whereas those due to tensor fluctuations decay away 
H. However the height of the Doppler peak also de- 
pends sensitively on several poorly determined cosmo- 
logical parameters (for example the Hubble constant and 
the baryon density) , on the ionization history of the uni- 
verse, and on theoretical model parameters (the details 
of the inflaton potential (pP]). These uncertainties make 
it difficult to unambiguously determine the size of the 
tensor contribution [pj. 

The polarization of the microwave sky could provide 
invaluable extra information, not far beyond the reach of 
current experiments. The degree of linear polarization 
expected from scalar and tensor fluctuations produced 
during inflation has been calculated ^-|l^, with the re- 
sult that for a given level of temperature anisotropy ten- 
sor perturbations do produce a somewhat larger degree 
of linear polarization. 

Recently we suggested using the temperature- 
polarization cross correlation (QT) as a further test [Q . 
Being linear in the polarization, this has some advantages 
for experiments in which noise in the polarization is lim- 
iting, for the latter averages to zero in {QT). The (QT) 
correlation also extends to larger angular scales than the 



auto-correlation function (QQ). It thus has particular 
relevance to experiments with large sky coverage, such 
as post-COBE satellites now being planned. Here we ex- 
tend our calculations of (QT) to the tensor case, with 
the result that a striking distinction emerges. On large 
angular scales the scalar and tensor (QT)'s have the op- 
posite sign, as a direct result of fundamental differences 
between scalar and tensor modes. Whether the rather 
small signal we predict is observable is unclear, depend- 
ing primarily upon the levels of foreground contamina- 
tion produced by our galaxy, which as we shall discuss 
below is still unknown. 

We begin with a simplified discussion. To first approx- 
imation, recombination may be treated as instantaneous 
and there is no polarization. The next approximation 
is to assume that each photon subsequently undergoes a 
single scattering, after travelling a comoving distance A. 
The photons we now receive from a given direction on 
the sky emanated from a sphere of radius A surrounding 
the scattering point (Fig. 1). The Thomson cross section 
is (Tt oc |e • e'l'^, with e and e' the initial and final photon 
polarizations. It follows that the photons we measure 
polarized along the y axis mostly come from above and 
below the scatterer, and those polarized along the x axis 
mostly come from the sides of the scatterer. 

In Figure 1 a single plane wave perturbation is shown, 
with its k vector perpendicular to the line of sight. For a 
scalar perturbation, the shading indicates the Newtonian 
potential <&, darker indicating more negative <&. Photons 
coming from the potential trough are redshifted, so it 
appears to us as a temperature trough on the sky. Pho- 
tons coming from above and below the scatterer fall into 
the potential trough before they scatter, and so are rel- 
atively blueshifted. Photons from either side suffer no 
such shift. So the net linear polarization is aligned par- 
allel to the temperature troughs. By rotating and su- 
perposing a set of such modes, one sees that tempera- 
ture cold spots are surrounded by a radial polarization 
pattern, hot spots by a tangential pattern. A simple 
quantitative estimate may be made in the long wave- 
length limit (fcA <IC 1), by averaging the Thomson cross 
section over the sphere of incident photons. The tem- 
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the radiation era. Using this, one finds Q 
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FIG. 1. A perturbation mode on the surface of last scat- 
tering, with the paths of photons which we see coming from 
a given point on the sky. For scalar perturbations, the di- 
rection of linear polarization (heavy lines) is aligned with the 
temperature troughs. For tensor perturbations the polariza- 
tion is instead perpendicular. 



perature anisotropy is ST/T « i<i>(s), with s denot- 
ing 'at the scatterer'. But the polarization depends on 
the quadrupole moment of $ on the sphere - if we ex- 
pand ^(x,y) in a Taylor series, the first terms that con- 
tribute involve second derivatives. One finds the frac- 
tional linear polarization measured using the x, y axes is 
Q = -^XHdl - dl){5T/T)[x,v), where T{x,y) is the 
two dimensional temperature field. 

Gravitational waves are traceless, transverse excita- 
tions of the metric, g^i, — a?'{ri^i, -\- h^^), where a is 
the cosmic scale factor. Consider a gravity wave with 
the same k vector and only h^z — ~hyy nonzero. The 
shading in Figure 1 indicates negative h^z- The tem- 
perature distortion is given by the Sachs- Wolfe integral 
~h J hij,on^n^ « ^hzzis), with the dominant contribu- 
tion coming from the path shared by all photons from 
the scatterer to us. So again Figure 1 shows a temper- 
ature trough on the sky. But for gravity waves initially 
outside the horizon, if hzz is negative, hzz is positive, so 
the proper distance is being expanded in the z direction, 
contracted in the y direction. Now we see a fundamen- 
tal difference with the scalar case. Photons arriving at 
the scatterer from above and below are unaffected (since 
the gravity wave is transverse), whereas photons coming 
from the two sides are blueshifted {hyy being negative). 
The induced polarization pattern is therefore perpendic- 
ular to the temperature troughs. Again by rotating and 
superposing one sees that for gravity waves the tempera- 
ture cold spots are surrounded by a tangential, hot spots 
by a radial polarization pattern. 

Repeating the 'single scattering' calculation explained 
above, we find Q « 20'^'''^^. o('^)- ^ gravity wave mode 
initially outside the horizon evolves as hzz oc 1 — Ak^rj^, 



dy){ST/T){x, y)Xriis, Vis being the conformal time at last 
scattering. 

If the temperature autocorrelation function is approx- 
imately scale invariant, then (T(6', (/))T(0)) - ln(l/6'), 
which is a rough approximation for both scalar and ten- 
sor cases for 6is <^ 9 <^ 1, 9is being the angle sub- 
tended by the horizon at last scattering. Using the re- 
lation between Q and T derived above, one sees that 
{Q{9,(f>)T{0)) ~ cos(2(/))0~^, with coefficients of opposite 
sign for the scalar and tensor cases. On smaller angular 
scales, one is sensitive to modes of wavelength smaller 
than the horizon, which have begun to oscillate by the 
time of last scattering. In the scalar case oscillations in 
the photon-baryon fluid density S cause the sign of (QT) 
to oscillate, since T (x S and Q on 6 are out of phase Q . 
Likewise in the tensor case, we have T oc h and Q oc h, 
so the sign of (QT) also reverses. However the redshift- 
ing away of gravity waves reduces these oscillations to 
negligible levels, so that {QT) is actually positive for all 
6* <C 1. At very small angles, from the cos(2(/)) depen- 
dence and analyticity it follows that (QT) vanishes like 
9\ 

To calculate (QT) accurately we evolve the photon dis- 
tribution function, f(x, p,t), using the relativistic Boltz- 
mann equation for radiative transfer with a Thomson 
source term, keeping terms to first order in the per- 
turbation. The distribution function, f, is a four di- 
mensional vector describing the intensity and polariza- 
tion degrees of freedom, with components related to the 
Stokes parameters: T,Q,U, and V ||l2|. The Boltzmann 
equation can be rewritten in terms of the perturbed dis- 
tribution functions, or brightness functions, defined as, 
A^ = 4(5/,/ (To^V where Tq is the mean GBR temper- 
ature, / is the unperturbed Planck distribution, Sfi is its 
first order perturbation, and i — T,Q,U or V. 

We evolve the coupled equations by expanding the per- 
turbations in in plane waves. In the case of scalar pertur- 
bations, one needs to evolve only two transfer equations, 
those for the components of f corresponding to Stokes' 
parameters T and Q g. In the tensor case the U com- 
ponent does not vanish but Polnarev [g| has shown with 
the proper choice of variables. 



A^(^, (/.k) = a(Ai)(f - m') cos(20k) 
A«(//, (/-k) = /3(Ai)(f + /i') cos(20k) 
A^(/i,0k)=/5(M)2Aisin(20k) 



(1) 
(2) 
(3) 



(where /i = k • q and (f>k is the polar angle of q about k) 
the four Boltzmann equations reduce to just two coupled 
equations: 



a + ikpia ~ h — aTnea[a + ^] 
/? -I- iki^P = -arUealP - *] 



(4) 
(5) 



where, 
* ^ ^ I d^,'[{l + 6m'^ + f,''}f3{f,') - (1 - /i'^)^a(/.')]. 

(6) 

Here, n^ is the density of free electrons. These equations 
are evolved, as in the scalar case 0], by expanding a 
and f3 in Legendre polynomials (i.e., a{fi) — X]i(2' + 
l)a;P;(/i)), converting them into a hierarchy of ordinary 
differential equations |g]. 

Once these variables are evolved to the present epoch, 
correlation functions are evaluated by summing over 
k and possible polarizations. One finds that the 
temperature-polarization cross correlation function is. 
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where {9, 4>) are the usual spherical polar angles of q, 
and the axes used to define the Stokes parameters are e^ 
and By. Here, A^y) = ij,'^a{n), B{iJ,) = (1 + fi'^)(3{fj,) 
and the constants aw and dw are given by aui = 
J_^dxPi{x)P^{x), and dw = J_^dxPi{x)P^{x) which 
have simple closed form expressions. 
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FIG. 2. QT correlation functions for tensor perturbations 
in a universe with standard recombination, and in a fully 
ionized universe, both normalized assuming that the entire 
COBE detection is due to gravity waves. Shown for com- 
parison is the result for scalar perturbations in a standard 
scenario [11]. 

Figure || shows (QT) for = 0. As for (QQ), on small 
scales the tensor (QT) is small in comparison to the scalar 
{QT). At large 9 , however, the signals have comparable 



magnitudes and opposite signs. With substantial early 
ionization, polarization is greatly enhanced on large an- 
gular scales (note also that a geometrical effect causes 
the tensor (QT) to reverse sign for 9 > 30°.) 




FIG. 3. 20° X 20° temperature map resulting from tensor 
perturbations, smoothed on 3°, with the correlated compo- 
nent of the polarization overlaid. 




FIG. 4. As above, but for adiabatic scalar perturbations. 

Since the primordial fluctuations are assumed to be de- 
scribed by gaussian statistics, the temperature and po- 
larization are completely described by (QT) (TT) and 
{QQ)-, the latter two calculated in p and [|0| respec- 
tively. Using these it is straightforward to construct real- 
izations of the microwave sky [[ll| . The total polarization 
Q is composed of parts which are correlated, Qc(q)i f^nd 
uncorrelated, (5c/(q), with the temperature anisotropy. 



Figure Q shows the correlated component overlaid on the 
temperature field. The length of each vector is propor- 
tional to [Q-^{d) + Uq{9)]2 and the orientation is given 
by 20 — tan~^{Uc/Qc)- Hot spots are seen to be asso- 
ciated with radial polarization patterns, while cold spots 
are surrounded by tangential polarization patterns. For 
scalar fluctuations on similar scales, however, the oppo- 
site is true (Fig. ||). 

The power spectrum of Q is compared to that for Qc in 
Figure 0, showing that in the tensor case the polarization 
is more strongly correlated with temperature than in the 
scalar case. The variance cr?, , proportional to the area 
beneath the curve, comprises more than one third of the 
variance of the total polarization, aQ. By comparison, 
for scalar fluctuations the correlated variance is barely a 
seventh of the total. 
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FIG. 5. The power spectra of the total (dashed line) and 
temperature-correlated (solid line) polarization in a universe 
with a standard thermal history normalized to the COBE 
measurement of the 10° smoothed temperature variance. 

Are these temperature-polarization correlations ob- 
servable? The main problem is whether sufhcient sky 
coverage can be obtained for a meaningful test. To give 
ourselves the best possible chance we consider an exper- 
iment which covers the whole sky, with very high signal 
to noise in the temperature measurement. A natural ob- 
servable is then C = ^ / dQQQc , being the integral of 
the 'predicted' polarization Qc computed from the tem- 
perature alone, multiplied by the 'observed' polarization 
Q. Since we are only able to make a single measurement, 
a prediction for O is only testable if S/N = O/ctq 3> 1. 
In Table I we have computed the signal to noise ratios for 
different beam sizes in a hypothetical full sky experiment. 
These numbers include only the effects of cosmic variance 
- instrument noise may be modeled by adding an ex- 
tra temperature-uncorrelated white-noise component to 
Q with variance a-j^ . This has the effect of reducing S/N 



by a factor (1 -f cr 
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which may be computed from 



the values for uq given in the Table. 

Apart from the forbidding challenge of building detec- 
tors with the [iK sensitivities required, the biggest ques- 
tion for using this technique to distinguish scalar and 
tensor perturbations is whether the polarization signal is 
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TABLE I. Total and correlated polarization in different 
theoretical scenarios for various beam sizes Q (FWHM). ctq 
and GQ^ are given in {iK, with the theories normalized to 
rms temperature fluctuations of 40/i-ft' for Q = 10°. 



swamped by foreground contamination from our galaxy. 
At frequencies below ~ 100 GHz, synchrotron emission 
is a significant background to CMB experiments, and 
a typical expectation is linear polarization at the level 
of 5 — 10%, significantly larger than the effects we are 
looking for (except perhaps in the fully reionized scenar- 
ios). At higher frequencies, dust emission becomes the 
main background of concern. Here the situation appears 
much more optimistic, with linear polarization around 
1% being typical ||l^ , although in special regions of very 
high magnetic field, it can be substantially higher p4[ . 
The frequency dependence of the dust emission may al- 
low one to subtract out dust-associated polarization us- 
ing multifrequency measurements. In any case we hope 
the prospect of observing the intriguing signals investi- 
gated here will stimulate further study of the likely back- 
grounds to future full sky polarization measurements. 
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